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The nonlinear dynamics of an obliquely oriented wave packet at sea surface is studied both 
analytically and numerically for various initial parameters of the packet, in connection with the 
problem of oceanic rogue waves. In the framework of Gaussian variational ansatz applied to the 
corresponding (1+2D) hyperbolic nonlinear Schrodinger equation, a simplified Lagrangian system 
of differential equations is derived, which determines the evolution of coefficients of the real and 
imaginary quadratic forms appearing in the Gaussian. This model provides a semi-quantitative 
description for the process of nonlinear spatio-temporal focusing, which is one of the most probable 
mechanisms of rogue wave formation in random wave fields. The system is integrated in quadratures, 
which fact allows us to understand qualitative differences between the linear and nonlinear regimes 
of the focusing of wave packet. Comparison of the Gaussian model predictions with results of direct 
numerical simulation of fully nonlinear long-crested water waves is carried out. 
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I. INTRODUCTION 

Anomalous waves at the ocean surface and in many 
other physical systems (known also as giant waves, rogue 
waves, freak waves) have been a popular subject of 
present-day scientific research (see, e.g., reviews 00, 
and references therein). The most frequently used math¬ 
ematical model for this phenomenon is the (1+1D) non¬ 
linear Schrodinger equation (NLSE) with the focusing 
nonlinearity. In particular, such equation describes the 
complex amplitude of the main harmonic of a quasi- 
monochromatic weakly nonlinear surface wave over plane 
potential flows of an ideal fluid Q. Due to nonlinear self- 
focusing, the modulation instability develops in a suffi¬ 
ciently lengthy and high wave group 00. As the re¬ 
sult, a one-dimensional anomalous wave arises. Such a 
scenario is confirmed by numerical and laboratory exper¬ 
iments The full integrability makes the (1+1D) 

NLSE attractive for analytical studies and provides many 
exact solutions which describe some important properties 
of real rogue waves (see, e.g., [T^Tti 1. 

In the Nature however sea waves are rather far from 
the one-dimensional model (see, e.g., fl7i-l23i]). When on 
the free fluid surface there are disturbances depending on 
both horizontal coordinates, the wave dynamics becomes 
more complicated already in the framework of the corre¬ 
sponding NLSE, not to speak about strongly nonlinear 
regimes. First, (1+2D) NLSE is not integrable. Second, 
in the case of deep-water gravity waves, the spatial dif¬ 
ferential operator in it appears hyperbolic: 

2iip t + %p xx - ip yy + \ip\ 2 ip = 0. (1) 
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Here the non-dimensional variables are in use: k$A* —»• 
ip, coot —> t, 2koX —> x, V^koy —> y, where A{x, y, t) is 
the complex envelope of the main harmonic, ko is the 
wave number, u>q = 2i t/T 0 = y/gko is the frequency of 
the carrier wave (the reference frame is moving with the 
group velocity u gr = (l/2)^/g/ko along a: axis). Although 
Eq.([l]) has some disadvantages, which are not discussed 
here, but it is more suitable for approximate analytical 
studies (because in all its modified variants, some pseudo¬ 
differential operators appear inevitably, difficult to calcu¬ 
late them “by hands”). 

Thus, in the transverse direction the nonlinearity acts 
in the defocusing manner, while in the longitudinal direc¬ 
tion it works for focusing of a wave packet. In addition, 
the linear dispersion makes its own defocusing contribu¬ 
tion. Owing to the indicated factors, the behavior of a 
nonlinear group of waves on the two-dimensional free sur¬ 
face of the fluid is more diverse as compared to the one¬ 
dimensional case. Depending on initial conditions, the 
nonlinearity can in some cases reinforce the linear focus¬ 
ing, but as well it can destroy the focusing in other cases. 
In the case of reinforcement, a rogue wave arises. Besides 
that, anomalous waves can appear as a result of interac¬ 
tion between previously formed coherent structures poj . 
However, for the usual oceanic wave fields where the pres¬ 
ence of coherent structures is hardly probable, more ac¬ 
tual seems the scenario when anomalous waves rise due 
to occasional spatio-temporal focusing (see jH, E5, and 
references therein). 

In general, the question about optimal conditions for 
the nonlinear focusing of water waves is still far from be¬ 
ing clear. It is only obvious that, similarly to the ID 
case, the initial wave group should be sufficiently high 
and/or extensive. But an effect of geometric shape of 
the packet, as well as phase modulation, on the develop¬ 
ment of anomalous wave is poorly studied. In the absence 
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of exact (1+2D) NLSE solutions, which could describe 
quantitatively all such processes, it seems reasonable to 
carry out approximate consideration of the dynamics of 
an idealized solitary wave packet characterized by a few 
time-dependent parameters. Thereby, the above qualita¬ 
tive reasons could be concretized, and a reference point 
could be put along the direction of development of future 
more accurate theory of three-dimensional rogue waves. 

The purpose of this work is to investigate the nonlin¬ 
ear dynamics of a wave packet in the framework of the 
full Gaussian variational ansatz, i.e. in the case when the 
quadratic form under the exponential has an off-diagonal 
part. Such a packet has an elliptical shape, with main 
axes oriented at some time-dependent angle with respect 
to the coordinate axes. Contrary to the more widely 
known diagonal ansatz (see, e.g., |26 MTm 1 ). the full Gaus¬ 
sian ansatz was previously applied to the elliptic NLSE 
only [35], [36). The present work fills this gap in the theory. 
Here a system of variational equations will be derived for 
parameters of the ellipse, and its full integrability will 
be demonstrated. The fact of integrability takes place, 
first, because the third power of the NLSE nonlinearity 
is accorded with the two spatial dimensions, and second, 
because Eq. (HD has a special integral of motion, namely 
the hyperbolic analog of the angular momentum, 

J — + - ipil>*)\dxdy = const. (2) 

The knowledge of general structure of analytical so¬ 
lutions of the variational model will allow us to under¬ 
stand qualitative differences between the linear and non¬ 
linear regimes of the wave packet focusing. Besides that, 
with the purpose of immediate testing of the Gaussian 
model, a comparison between its predictions and the re¬ 
sults of numerical simulations of fully nonlinear long- 
crested waves will be carried out. We shall see that the 
agreement can be rather good even for strongly nonlinear 
regimes, when the very NLSE is not applicable already. 


II. VARIATIONAL ANSATZ 


Let us first remind that in the simplest — diagonal 
variant, the Gaussian ansatz describes an elliptical 
wave packet having the symmetry axes coinciding with 
the coordinate axes (see, e.g., |26t - l3l| l: 
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( 3 ) 


Besides the longitudinal size X(t) and transverse size 
Y (t ), such a packet is characterized by four more real 
quantities: U, V, N, and 0. Parameter N does not de¬ 
pend on time because 4irN = f \ijj\ 2 dxdy is an exact inte¬ 
gral of motion of NLSE, while the conjugate variable 0 is 
cyclic. Parameters U and V describe a phase modulation, 
with positive/negative values of U or V corresponding to 
defocusing/focusing configurations respectively in x or y 


direction. Here it is necessary to note that for applica¬ 
bility of NLSE, the condition of narrow spectral width of 
the wave packet should be satisfied. In our case it prac¬ 
tically means the following: (X,Y) > 10, (U,V) < 0.1. 
Equations of motion determining the temporal behavior 
of the unknown quantities X, Y, U, and V, are derived 
by substitution of the variational ansatz © into the La- 
grangian functional of the NLSE: 

C = J\itptV ~ iip^t - \4>x\ 2 + \4> y \ 2 + 10| 4 /2 )dxdy. (4) 


In this way we obtain a reduced Lagrangian C depending 
on TV, X(t), Y(t), U(t ), and V(t). The action integral 
J Cdt should be then variated. The main outcome of 
this standard procedure is the following. Variation of the 
reduced action gives in particular that U = X. V = Y. 
Finally we arrive at the following system of differential 
equations of the Newtonian type: 

-IN 1 TV 

X= X*~X 2 Y' V= Y^ + Y 2 X' (5) 

Note that equations © constitute a Lagrangian system 
with the Lagrangian functional L oc C/N , where 


2 L = X 2 ~Y 2 - 


X 2 + Y 2 


2TV 

XY’ 


( 6 ) 


so the kinetic energy is not positively-definite (the second 
“particle” has a negative mass). 

It should be said that equations of the type © and 
their generalizations are actively used in plasma physics, 
in nonlinear optics, and in atomic physics to investigate 
behavior of NLSE and Gross-Pitaevskii equation solu¬ 
tions in two and three spatial dimensions (see papers 
26-134]. and references therein). In particular in work 
271 ]. analytical and some numerical solutions were pre¬ 
sented namely for the hyperbolic (1+2D) NLSE, as in our 
case, but with such parameter values that are typical for 
nonlinear optics, not for water waves. In relation to the 
problem of anomalous sea wave focusing, this simplified 
variational ansatz was applied only very recently in the 
author’s work [37]. 

Let us now turn our attention to the full Gaussian 
ansatz, i.e. when off-diagonal elements of the quadratic 
forms are involved. For further consideration, it is con¬ 
venient to introduce linear combinations of the spatial 
coordinates, corresponding to a hyperbolic rotation with 
the parameter y (analogously to the elliptic NLSE case, 
but with the difference that there the rotation was usual 
trigonometric): 

x = a; cosh % + j/sinhx, y = x sinh x + y cosh(7) 

Instead of expression © we now write 


4TV 

V = \l Yy 6XP 


y 2 Ux 2 Vy 2 

-+i —i —i^xy+icp 


2X 2 2Y 2 2X 
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( 8 ) 
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Strictly speaking, with \ ^ 0 the parameters X and 
Y are no longer the longitudinal and transverse sizes 
of the packet, but we shall continue to call them so for 
brevity. Parameter \ cannot be too large, otherwise the 
wave packet becomes spectrally wide and goes beyond 
the NLSE applicability domain. Practically, as it will be 
seen from the results of direct numerical simulation of 
nonlinear waves with realistic parameters, the inequality 
|x| < 0.5 should take place. 

Note that Gaussian form © is one of exact solutions 
of the corresponding linear Schrodinger equation, with 
appropriate temporal dependencies of the parameters ap¬ 
pearing there. Consequently, at small N the approxima¬ 
tion © is certainly valid. In practice, values N ss 2-4 are 
of interest. In this case, generally speaking, only quali¬ 
tative agreement of the results of variational model © 
with the solutions of the NLSE (and, the more so, with 
the completely nonlinear dynamics of waves on water) 
can be expected. In more detail the comparison between 
the variational and numerical solutions will be discussed 
some later. An undoubted benefit of this approxima¬ 
tion is that it provides a semiquantitative description of 
the process of spatio-temporal focusing, which is one of 
the most probable mechanisms of the formation of rogue 
waves under real conditions (see [24], HH, and references 
therein). 

Having substituted the trial function © and its partial 
derivatives (expressed it terms of x and y) into the La- 
grangian after straightforward calculations we obtain 
that variable \ is cyclic, while the corresponding con¬ 
served quantity is (taking into account that N = const) 

M = 7 ( X 2 + Y 2 ) = const. (9) 

Herewith, the temporal dependence of parameter \ is 
determined by the equation 

X = M(Y 2 - X 2 )/(X 2 + Y 2 ) 2 . (10) 


III. INTEGRABILITY OF THE MODEL 

An important advantage of system m is an exact 
integrability, because here in analogy with the elliptic 
case, the separation of variables turns out to be pos¬ 
sible. Integrability of the variational approximation is 
typical namely for NLSE in 2D space, because in 3D the 
reduced Lagrangian instead of 2 N/ (XY) would contain 
2TV/ (XYZ ), and it would destroy the homogeneity of the 
equations. It is well known from the analytical mechan¬ 
ics that for integration of the above system one should 
use the hyperbolic coordinates: 

X = Q sinh £, Y = Q cosh£. (14) 

These coordinates describe sector X < Y which is most 
interesting in application to the problem of oceanic rogue 
waves, because a length of their crest usually 3-10 times 
exceeds the wave length Ao = 27r/fco, while the longitu¬ 
dinal packet size in the minimum reaches approximately 
one Ao (that corresponds to dimensionless values X m i n ~ 
6-10). The Lagrangian (fill) takes the form 

2 L M = Q 2 i 2 -Q 2 -^, (15) 

where function F(£) is defined by the formula 

4 4N M2 

sinh 2 (2^) sinh(2f) cosh 2 (2^) 

The corresponding equations of motion are 

Q + Qi 2 + ^r = °- ( 17 ) 

• • ■■ F'(f) 

2 QQi + Q 2 i + ^0- = 0. (18) 


The conserved quantity M is a hyperbolic analog of the 
angular momentum. As previously, we have U = X, 
V = Y, but the dynamics of unknown functions X(t) 
and Y (t) is described now by the following Lagrangian: 


2 L m = X 2 —Y 2 — 


X 2 Y 2 


2 N M 2 (X 2 -Y 2 ) 
XY + (X 2 +Y 2 ) 2 ' 


This Lagrangian is a hyperbolic analo g of the elliptic vari¬ 
ational problem considered in works ©I ©J. Equations 
of motion that follow from here, have the form 



N 

X 2 Y 


M 2 X 

(X 2 + Y 2 ) 2 


1 


o (X 2 — Y 2 ) 
2 (X 2 + Y 2 ) ’ 


( 12 ) 


After multiplying Eo. flTSl) by 2Q 2 £, we obtain in its left 
hand side a total time-derivative, so the integral of mo¬ 
tion follows: 


Q 4 f + F(C) = C = const. (19) 

Looking at the expression (HU), it is easy to understand 
that at sufficiently small negative values of the integra¬ 
tion constant C , equation F(£) = C determining “turn 
points”, has two real positive roots (£i and £ 2 ), with 
the motion of variable £ taking place between them. At 
C > 0 the “turn point” is sole. 

Eq. fTTl) is simplified to Q — —C/Q 3 . Its general solu¬ 
tion is written in explicit form as 



N 

Y^X 


M 2 Y o (Y 2 — X 2 ) 

(X 2 + Y 2 ) 2 ~ 2 (X 2 + Y 2 ) 


(13) 


The finding solutions of this system and the comparison 
of the corresponding evolution of wave packet with a di¬ 
rect numerical modeling of nonlinear waves is the essence 
of the present work. 


Q(t) = y/l(t-u) 2 -c/l, ( 20 ) 

where f* and I are two more integration constants. It 
is easy to check that I = —2 E, where E is the energy 
integral of the system: 

E=(Q 2 i 2 ~Q 2 + F(0/Q 2 )/ 2 . ( 21 ) 
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With taking into account Eq. (l20l) . it follows from Ea. (fl9l) 
that 


f _ d l _ = ± _ ( 22 ) 

where £o is the fourth — the last - integration con¬ 
stant. Formulas (1201) and (l22ll provide the full solution 
of the problem. In particular, with I > 0, C < 0, the 
dependence (1^1) can be represented in the following form, 


sinh 2£ (t) = P (ao + 


1 

, arctan 


i \)’ 


(23) 


where an even a-periodic function P(a) (having a min¬ 
imum at a = 0) is the inverse with respect to elliptic 
integral (here the integration variable z = sinh 2^): 


P ( ± f 

J S i 


zdz 


2^/(1 + z 2 )(Cz 2 + 4Nz - 4) - M 2 z 2 


= s, 


(24) 

and Si is the smaller root of the equation (1 + s 2 )(Cs 2 + 
4 Ns — 4) — M 2 s 2 = 0. The amplitude of wave packet 
equals to \JAN/XY, and it is given by the formula 


A(t) = V8N/^[I(t - f*) 2 - C/I\ sinh2£(t). (25) 


IV. SELECTION OF THE SOLUTIONS 



t/T 0 


FIG. 1: The temporal dependencies X(t) and Y(t) for N = 
3.0, M = 0, V 0 = 40, To = 50, U 0 = -0.05, Vo = -0.10. 
The solutions for a linear packet with the same initial condi¬ 
tions are also shown for comparison. What is important for 
formation of a rogue wave, the nonlinear focusing along the 
longitudinal direction goes here with acceleration for most of 
the time, while in the linear case it occurs with deceleration 
for all the time. The value of X in the minimum has been too 
small for validity of the NLSE, but they are initial conditions 
of this type which in reality result in formation of anomalous 
waves of the limiting steepness. 


In fact, we shall not use the analytical solutions, ex¬ 
pressed in terms of special functions, for finding the de¬ 
pendencies X(t) and Y(t) with some given initial values. 
Practically it is much simpler and faster to carry out 
highly accurate numerical simulation of the variational 
equations m and m immediately. One should also 
remember that the absolute accuracy in solving approxi¬ 
mate variational problem is hardly necessary. The knowl¬ 
edge of the analytical structure of the solutions will be 
needed for a general understanding of properties of the 
dynamical system under consideration. 

We are interested mainly in those solutions from the 
whole diversity (BUI) and (E2l) . which result in essentially 
higher wave amplitude as compared to the maximally 
possible linear value A max im = y / 4./Vt/_ 00 V r _ 00 , where 
[/_oo and V-oo are the asymptotic values of the focus¬ 
ing parameters at large negative times. In essence, only 
such events are deserving the terrible title “rogue waves”. 
What is important, the most interesting solutions corre¬ 
spond to values of the parameters C « 0, M ~ 0. The 
condition C = 0 means that Q 2 = Y 2 — X 2 = I(t — t*) 2 , 
i. e. the plots of functions X(t) and Y(t) touch each 
other at t = f*. It is easy to make sure of the above 
assertion by numerical simulations of Eas. (fl2l) - (fT3l) [with 
fixed N, Uq, Vq, Y 0 , but different and M], and sub¬ 
sequent comparison of the maximal amplitudes. Typical 
dependencies X(t) and Y(t) for such a case are shown 
in Fig[T| Let us also note that the minimal value of the 
ratio (V/F) m i n = tanh£i, which is reached closely to the 


moment of maximal amplitude, at M = 0 and C = 0 can 
be found from the simple condition sinh 2£i = 1/iV, as it 
is clear in view of formula (1241) . 

It is worth noting that at M = 0 the cyclic variable 
X remains constant, as it follows from Eq. m- As we 
shall see soon, the physical wave fields corresponding to 
different values of y, strongly differ between each other. 
In particular, at y = 0 the wave picture is symmetric 
with respect to sign change of y coordinate, while at 
y / 0 there is no any symmetry (see FigsHE]). But 
in the framework of NLSE, the envelope of one solution 
coincides with the envelope of the other solution after 
a hyperbolic rotation of the dimensionless coordinates. 
The presence of this approximate symmetry significantly 
clarifies the question about optimal conditions for the 
formation of anomalous waves. 


V. COMPARISON WITH NUMERICAL 
EXPERIMENTS 

It is clear that any simplified model without verifica¬ 
tion remains just a toy. In order to make the status of the 
variational approximation more distinct, it is necessary 
to compare it with a more accurate model. In our case, 
numerical experiments were carried out in the framework 
of the fully nonlinear long-crested wave model which is 
described in detail in Refs.[38l [39|. The computational 
domain was a square of the size 5 km, with the periodic 
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a) 1=00 min 00.0 s. Surface elevation [m] 
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b) 1=11 min 15.5 s. Surface elevation [m] 



1 1.5 2 2.5 3 3.5 4 
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b) 1=11 min 15.5 s. Surface elevation [m] 
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FIG. 2: Wave picture at x — 0.0 (the other parameters 
are the same as in Fig.l): a) wave packet at the initial time 
moment; b) at the moment when the anomalous wave has 
risen. The crest length of the high wave here has been much 
longer than it would be in the linear case, which observation is 
in accordance with the Gaussian model prediction. Deviations 
from the Gaussianity are clearly seen, basically — in the form 
of a small tail consisting of longer waves. 

boundary conditions. The length of the carrier wave was 
chosen Ao = 125 m. The simulations were carried out for 
different sets of initial parameters of the packet. 

In Fig.2a and Fig.3a, examples are shown of how Gaus¬ 
sian initial conditions look. The “portraits” of anomalous 
waves developing from those conditions, are presented 
in Fig.2b and Fig.3b. In general, deviations from the 
Gaussian shape are clearly seen there. Despite the non- 


FIG. 3: Wave picture at x — 0.3 (the other parameters 
are the same as in Fig.l): a) wave packet at the initial time 
moment; b) at the moment of the anomalous wave rise. De¬ 
viations from the Gaussianity are even more prominent here 
than at x = 0.0. 


Gaussianity of the arising rogue waves, Figs.4-5 demon¬ 
strate that there is a qualitative agreement between the 
variational theory and the numerical experiment in the 
important aspect as the temporal dependence of the max¬ 
imal amplitude of the wave packet. It is seen that Gaus¬ 
sian model somewhat overestimates the height of anoma¬ 
lous wave (especially as parameter x increases to values 
about unity, when the approximate symmetry of the hy¬ 
perbolic rotation becomes invalid), and it moves to the 
future the moment of the maximal rise. However, it is 
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t/T 0 



t/T 0 


FIG. 4: Comparison of the temporal dependencies of the max¬ 
imal amplitudes for wave packets with different the other 
parameters being the same as in Fig.l. In each figure shown 
are numerically found the height of the highest crest and the 
depth of the deepest trough, and also predictions of the Gaus¬ 
sian models - linear and nonlinear. The numerical curves 
have an oscillating character, which property is related to the 
approximately two-fold difference between the phase veloc¬ 
ity of the carrier wave and the group velocity of the packet 
motion. Heights of crests are typically larger than depths of 
troughs, due to the presence of higher harmonics. Therefore 
the quantity that should be compared to the theoretical pre¬ 
dictions is a half-sum of envelopes of the numerical curves. 




I/Tq 


FIG. 5: Comparison of the maximal amplitudes for different 
sets of initial data than in Fig.4. It is seen that with increase 
of N, the time delay of the predicted peak in comparison with 
the numerical simulation increases. 


hardly reasonable to expect a detailed accordance from 
the simplest variational approximation. In any case, the 
time period of existence of anomalous wave is predicted 
by the Gaussian theory rather well. Also another im¬ 
portant effect is confirmed, the increase of the transverse 
packet size as compared to the linear theory (see Fig.l). 
This effect contributes to the subsequent more fast de¬ 
crease of the amplitude in comparison with the linear 
case. 


VI. CONCLUSIONS 

Thus, based on the variational analysis and on the 
comparison with the results of direct numerical simu¬ 
lations, it is natural to put forward the following the¬ 
sis. In the main approximation, a typical long-crested 
three-dimensional oceanic rogue wave at the nonlinear 
stage is, from the mathematical point of view, a dynam¬ 
ical system with just two degrees of freedom: X(t) and 
Y(t). This system contains two constant parameters (N 
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and M), which are related to the integrals of motion of 
1+2D NLSE and characterize a given group of waves. 
The system is integrable. Real sea states where one can 
expect anomalous waves, require values N ss 2-4. The 
most favorable focusing parameters for the rogue wave 
formation are: 1) M « 0; 2) on the initial stage of (oc¬ 


casional) focusing at some time moment t* the relations 
X * sa Y* ~ 20 — 40, —X* « — Y* > 0.05 should take place. 
To estimate probability of realization of such conditions 
in random weakly nonlinear wave fields is the task on the 
future. 


[11 C. Kharif and E. Pelinovsky, Eur. J. Mech. B/Fluids 22, 
603 (2003). 

[2] K. Dysthe, H. E. Krogstad, and P. Muller, Annu. Rev. 
Fluid Mech. 40, 287 (2008). 

[3] M. Onorato, S. Residori, U. Bortolozzo, A. Montina, and 
F. T. Arecchi, Physics Reports 528, 47 (2013). 

[4] V. E. Zakharov, J. Appl. Mech. Tech. Phys. 9,190 (1968). 
[51 T. B. Benjamin and J. E. Feir, J. Fluid Mech. 27, 417 

(1967). 

[6] V.E. Zakharov, A. I. Dyachenko, and O.A. Vasilyev, 
Eur. J. Mech. B/Fluids 21, 283 (2002). 

[7] A.I. Dyachenko and V.E. Zakharov, .JETP Letters 81, 
255 (2005). 

[8] V.E. Zakharov, A. I. Dyachenko, and O. A. Prokofiev, 
Eur. J. Mech. B/Fluids 25, 677 (2006). 

[9] V. P. Ruban, JETP Letters 95, 486 (2012). 

[10] A. Chabchoub, N. P. Hoffmann, and N. Akhmediev Phys. 
Rev. Lett. 106, 204502 (2011). 

[11] A. Chabchoub, N. Hoffmann, M. Onorato, and N. 
Akhmediev, Phys. Rev. X 2, 011015 (2012). 

[12] N. Akhmediev, V. M. Eleonskii, and N.E. Kulagin, Sov. 
Phys. JETP 62, 894 (1985). 

[13] N. Akhmediev and V. Korneev, Theor. Math. Phys. 69, 
1089 (1986). 

[14] N. N. Akhmediev, V. M. Eleonskii, and N. E. Kulagin, 
Theor. Math. Phys. 72, 809 (1987). 

[15] N. Akhmediev, A. Ankiewicz, and J. M. Soto-Crespo 
Phys. Rev. E 80, 026601 (2009) 

[16] M. Erkintalo, K. Hammani, B. Kibler, C. Finot et al., 
Phys. Rev. Lett. 107, 253901 (2011). 

[17] M. Onorato, A. R. Osborne, and M. Serio, Phys. Rev. 
Lett. 96, 014503 (2006). 

[18] M. Onorato, T. Waseda, A. Toffoli, L. Cavaleri et al., 
Phys. Rev. Lett. 102, 114502 (2009). 

[19] V.P. Ruban, Phys. Rev. E 74, 036305 (2006). 


[20] V.P. Ruban, Phys. Rev. Lett. 99, 044502 (2007). 

[21] D. H. Peregrine, Adv. Appl. Mech. 16, 9 (1976). 

[22] I. V. Lavrenov and A. V. Porubov, Eur. J. Mech. 
B/Fluids 25, 574 (2006). 

[23] V. I. Shrira and A. V. Slunyaev, Phys. Rev. E 89, 
041002(R) (2014). 

[24] C. Fochesato, S. Grilli, and F. Dias, Wave Motion 44, 
395 (2007). 

[25] V. P. Ruban, JETP Letters 97, 686 (2013). 

[26] M. Desaix, D. Anderson, and M. Lisak, J. Opt. Soc. Am. 
B 8, 2082 (1991). 

[27] L. Berge, Phys. Lett. A 189, 290 (1994). 

[28] L. Berge and J. J. Rasmussen, Phys. Rev. E 53, 4476 
(1996). 

[29] L. Berge, J. J. Rasmussen, E. A. Kuznetsov, E. G. 
Shapiro, and S. K. Turitsyn, J. Opt. Soc. Am. B 13, 
1879 (1996). 

[30] L. Berge, Phys. Reports 303, 259 (1998). 

[31] V. M. Perez-Garcia, H. Michinel, J. I. Cirac, M. Lewen- 
stein, and P. Zoller, Pliys. Rev. Lett. 77, 5320 (1996); 
Phys. Rev. A 56, 1424 (1997). 

[32] Yu. Kagan, E. L. Surkov, and G. V. Shlyapnikov, Phys. 
Rev. A 55, R18 (1997). 

[33] Y. Castin and R. Dum, Phys. Rev. Lett. 77, 5315 (1996). 

[34] F. Haas, Phys. Rev. A 65, 033603 (2002). 

[35] A. S. Desyatnikov, D. Buccoliero, M. R. Dennis, and Yu. 
S. Kivshar, Pliys. Rev. Lett. 104, 053902 (2010). 

[36] J. Abdullaev, A. S. Desyatnikov, and E. A. Ostrovskaya, 
J. Opt. 13, 064023 (2011). 

[37] V. P. Ruban, JETP Letters 100, 751 (2014). 

[38] V. P. Ruban and J. Dreher, Phys. Rev. E 72, 066303 
(2005). 

[39] V. P. Ruban, Eur. Phys. J. Special Topics 185, 17 (2010). 



